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§ ! Abstract 
CN 

The decoherence effect due to emission of gravitons is examined. It shows the 
same qualitative features of the QED effect which has already been investigated, it is 
obviously much weaker, wholly universal and shows a stronger energy dependence. 
The result can be extended to photons, they also may undergo decoherence due 
to graviton emission. For this limited aim the incomplete status of the quantum 
gravity, in comparison with QED, is not source of severe difficulties because all the 
O- ■ effects are attributed to the infrared sector of the dynamics. 
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I. Introduction 

The loss of quantum coherence due to bremsstrahlung process has been stud- 
ied in detail [1,2] for the electromagnetic case. The electromagnetic radiation is 
certainly present whenever charged particles undergo scattering processes, so, with- 
out assuming that the radiation is the unique process giving rise to the loss of 
quantum coherence, we may be sure that electrodynamical effects are there. It is 
also evident that this kind of decoherence mechanism makes a sharp distinction 
between charged and uncharged particles, it becomes therefore natural to look for 
other kinds of radiation that may implement the same effect in a universal way. 
The obvious candidate is the gravitational radiation since the gravity is the most 
universal form of interaction known at present, so in this note the effects of gravi- 
tational bremsstrahlung are examined. Some possible source of difficulties, appears 
immediately: when we investigate the effects of photon emission we have at our 
disposal QED which probably the best settled (and the oldest) branch of quantum 
field theory, the quantum theory of gravitation still waits a final shape. The small 
sector of the general quantum mechanics of gravitation that concerns the infrared 
radiation can be treated in close analogy with QED as it can be seen from the gen- 
eral treatment of the infrared radiations given by Weinberg [3] . This treatment will 
be systematically followed in this note, some points will be recalled or re-elaborated 
when necessary. The difference that cannot be neglected is the fact that the emis- 
sion of photon is due to the electric charge which is a Lorentz scalar and that it 
is quite possible to consider situations where every interacting particle keeps its 
charge, the emission of gravitons is due to the energy-momentum four vector and 
that this quantity is certainly not kept by the single interacting particle although it 
is conserved in the overall process. This fact is conceptually obvious but gives rise to 
some technical difficulties: for this reason in this paper the perturbative treatment 
is presented formerly, in sect. II because there the effects of the overall conservation 
are better shown, after an extension of the Bloch-Nordsieck model in proposed in 
sect. Ill where one sees that the effect of decoherence can be studied following each 
particle individually. The fact that also photons may undergo decoherence due to 
gravitational bremsstrahlung is noted. 

II. Infrared gravitons 

A. Infrared compensation for standard states 

Some results given in ref. [3] will be explicitly needed, so they are now recalled 
and adapted to the particular problem under investigation. Since the source of the 
gravitons is the energy-momentum tensor one must consider the whole scattering 
process, where this source is conserved, it is not enough to follow a particular line as 
it is allowed in QED [5]. To avoid unessential complications the scattering process 
is an elastic two-body collision, the masses are taken to be equal. 

The units c = h = 1 and the coupling k = V8nG are used (G is the Newtonian 



2 



constant), the Minkowski metric is denoted by rj^, the Mandelstam invariant s, t, u 
are also employed. 

If one calls M the amplitude for the scattering without further radiation 
Fig.l 

then one obtains the amplitude for the scattering with the radiation of a soft 
graviton [3]: 
Fig.2 



M r = M ■ i« 



q ■ k q' ■ k p ■ k p' ■ k 



{^)-" /2 f,u (1) 



One can get this expression working out e.g. the soft limit of the emission by a 
spinorial particle, using the vertex [6]*: 

\iK{\[^{ Pl + P2 )» + Y(Pi +V^[(P1 +P2) - 2m]} 

but it is quite general, it does not depend on the spin of the emitting particle. The 
polarization vector is and since the source is conserved the sum over polarization 
is simply given by: 

^ ^ f puf per ^-pu per 

i 

So the transition probability, with emission of an infrared graviton is given by 



i .. , ., r d s k „ x 4 / i i i i x 

'' ' \(q-k) 2 (q'-k) 2 (p-k) 2 (p'-k) 2 ) 



2(q ■ q') 2 — m 4 ^ 2(p ■ p') 2 — m 4 2(q ■ p) 2 — m 4 ^ 

q ■ kq' ■ k p ■ kp' ■ k q ■ kp ■ k 

2(q-p') 2 -m 4 2(p • q') 2 - m 4 2(p' ■ q') 2 - m 4 

q ■ kp' • k p ■ kq' • k p' • kq' • k 

The notation X\ will be used to indicate the sum of the first four terms, while X 2 will 
be used to indicate the sum of the last six terms. The integration over the energy 
of the radiated graviton goes from a minimum A, and a maximum A that has no 
role in the discussion, but for allowing the low-energy approximation. Performing 
the integration the following result is obtained: 

' \M \ 2 K 2 m 2 [l + D{p-p/m 2 )-D{p-q/m 2 )-D{p'-q/m 2 )\\n{K/\) (4) 



2tv 2 
where 



2x 2 — 1 

D{x) = arccoshx D(l) = 1 D(x) — > 2x \rv2x for x — > oo. (5) 

Vx 2 — 1 



In that paper the k is twice the k used here. 



3 



Now one must look at the virtual corrections: there are six terms corresponding to 
the correction in the channels s,t, u respectively, 
Fig.3 

working always in the i.r. limit the corresponding amplitude may be expressed 

as: 

M v = 2[M S + Mt+M u ] (6) 
and the correction, at the order k 2 is given by 2^t{M M.%) 

Ms = t^i^Mo I -f— t v " \\, 2 ''\ . , (7a) 
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Mt = 7S34«^o / TT-- ~ , ~ (76) 



(2tt)^ 
i 

• M » = <2^ k2a< ° y ^%f^i <7c) 

The integration over the virtual momenta which takes into account the contribu- 
tion of the photon pole gives the result needed to cancel one part of the real i.r. 
divergence [4,5], the term X2 ■ The contribution of the second pole of the s-channel 
term is easily calculated in the frame $ = 0, where it corresponds to l a = — ie with 
the result 

M f> == ^m4(p ■ P 'f - 1". 4 ]^ J j 1 = in m 

This term is imaginary, so it does not contribute to the correction to order k 2 , it 
corresponds (in QED) to the perturbative expansion of the Coulomb phase [4, 7], here 
it correspond to the gravitational elastic rescattering of the incoming or outgoing 
particles. 

There are still some other contributions, that are needed to cancel the divergent 
terms in X\. Here a difference with respect to QED is found , in fact in QED either 
one chooses the Yennie gauge [8] , so that the virtual corrections till now considered 
will completely cancel the i.r. divergent part coming from real radiation, or, when 
using the Lorentz gauge, one look at the i.r. divergent part of the renormalization 
constant Z 2 for the wave function of the charged particle. Here there is no sim- 
ple analogous of the Yennie gauge and a general renormalization procedure is not 
available. However when a self energy correction is inserted with the procedure: 

1 -sCp)^ 1 



p — m p — m p — m 

one must keep in the limit p — > p p 2 = m 2 , i.e. when the insertion is performed 
on one external leg, the correct position of the mass pole and the correct residuum 
at the pole, 
Fig.4 
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this amount to two subtractions that, in the i.r. regime are well defined; in 
particular the second subtraction term contains the factor <9£/<9p|^^_, when the 
derivative acts on the numerator of E it does not yield i.r. divergences, because £ 
itself is not i.r. divergent, on the denominator the derivative acts in the same way as 
in QED, and yields therefore the divergent terms needed to cancel the term X\. The 
important, and unfortunate, difference is that now these and similar subtractions 
does not avoid the ultraviolet divergences, but for our limited aim only the infrared 
sector matters. 

In simple formulation the result is that the i.r. radiation has a matrix element 

M r (k) = M ■ K(3(k) (8) 
The virtual correction has an i.r. contribution 

M v = M ■ n 2 p v (9) 

and that 



/ 



rl 3 k 

\[3(k)\ 2 ^ + ^Pv = C (10) 



The term C is finite, in the true i.r. limit is zero. 
B. Superposition of two states of motion 

If the two-particle final state is the superposition of two states of motion, like 
\T >= [\qi,q[ > +|</2, Q2 >]/v^2, then the transition probability without any soft 
particle correction is: 

±\mP + mP\ 2 (11) 

The transition probability with emission of one soft particle of momentum k is: 

\ \M Q X) ■ (k) +M a 2) - (k)\ 2 (12a) 
The transition probability including the correction for one virtual soft particle is: 

\\M^ ■ (1 + K 2 p^) + M 2) ' (1 + ^ 2 P[ 2) )| 2 (13a) 
so the virtual correction to the order k 2 may be written as: 

C v = K'QAtfW) + lM 2) | 2 p[ 2) ) +2K(M 1) M 2) *)(P1 1) +Pi 2) ) (136) 
while the total infrared real correction is: 

(126) 
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By using the compensation discussed previously between real and virtual cor- 
rection we find that the only a surviving addendum contains the interference term. 

C = C r + C v = -K 2 M(MPM*W) j \jp\k) - fP\k)\ 2 ^ (14) 
This expression can be recast in a form similar to eq.(4) 

C = -2nK 2 m 2 ^(M^M:^) 

[1 + D( qi ■ q'Jm 2 ) - D( qi ■ q 2 /m 2 ) - D( qi ■ q' 2 /m 2 )] ln(A/A) . 

If the two final states \qi,q[ > and 192,92 > which build up the superposition are 
not very different, it is convenient to use the expansion: 

q2 = qi + Sq q' 2 = q[ - 5q 

In the centre-of-mass frame the common three momentum is Q 2 = \s — m 2 and 
Sq ■ qi = —Sq ■ q[ = 2Q 2 sin 2 0/2, being the angle, that is now assumed small, 
between the two final directions q 2 and q± 

C = -A-kA{M^ M* {2) )Q 2 sin 2 0/2 [D(p ■ p'/m 2 ) - D(l)] ln(A/A) . (16a) 

From previous expressions 

£(1) = § D(x) -> 2(1 + In 2x) for x -> oo (166) 

When the radiating particle is massless some minor modification are required: if 
m — > then x — > oo in eq. (5) so the asymptotic form of D(x) is certainly correct 

C = -2tvk 2 ^(M i^Ml^) [m 2 + 2 Ql ■ q[ ln(2 qi ■ q[/m 2 )- 

2q x ■ q 2 ln(2gi • q 2 /m 2 ) - 2 Ql ■ q' 2 ln(2 qi ■ q 2 /m 2 )] ln(A/A) 



or also 

C = - ATtK 2 ^{M^M* ^)([ qi ■ q[ ln(2 qi ■ q'Js) - Ql ■ q 2 ln(2 gi • q 2 /s)- 
Qi ■ q' 2 H 2 Qi ■ v'2/s)] + [(Qi • 9i - 9i • 92 - 9i • q'2) ln(s/m 2 )]) ln(A/A) 



(17) 



Energy-momentum conservation gives gi + q[ — q 2 — q' 2 = so that in the limit 
qf = m 2 — > the second parenthesis goes to zero. 

In analogy with eq. (16) it is possible to give the expression for small 

C = 4ttk 2 3?(A^ 1) .M* (2) )2Q 2 sin 2 0/2 [2 In sin 0/2 - 1] (18) 

Comparing this expression with eq.(16a) it appears that the dependence on the 
presence of a logarithmic term is the only remnant of the massless condition, at least 
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for small angles. The fact the it is possible to deal with the i.r. divergence coming 
from massless particles emitted by a massless source in a simple way is, as shown in 
[3] , a peculiarity of the graviton emission, with photon emission the situation would 
be much worse; experimentally we know that there are massless particle having 
gravitational interactions, the photons, but there are not massless charged particles 
and if we look for charged massless particles in quantum chromodynamics we find 
indeed a very complicated infrared behaviour. 



III. Explicit description of the time evolution 
A. The Bloch and Nordsieck model reviewed 



The usual tools of scattering theory with emission of soft massless bosons going 
back to Bloch and Nordsieck model will be used [4], in that form which has been 
reviewed in the previous treatment of the QED effects[2], so only the points that 
show significant differences will be presented. The main difference is that, as dis- 
cussed in section II, it is not obvious that it is enough to follow only one particle 
in its fly, if we look to the final state, one should take into account both outcoming 
particles and together with their radiation and see how the compensation between 
real and virtual corrections is realized in this case. The following analysis shows, 
however, that the compensation happens in transparent way so that the consider- 
ation of only one particle, together with its radiation is, at the end, justified. The 
scattering process (for particles of equal masses) is described by the Hamiltonian 
[3]: 

2 

H = H Q + V = H + J2 V J + A ( 19a ) 

j=i 

H = -ivi ■ di + — - iv 2 ■ d 2 + — + / cPkual a k „ (196) 
7i 72 ^ J 

V> = -(^p7a£ / ^74^V^-^ • (We) 

The polarization tensors for the gravitons are expressed here in non-covariant Coulomb 
gauge and satisfy the well-known relations 

^ ] ffcjyffcjy = \\_ u ac u bd + u ad u bc ~ u ab u cd] u ij = $ij ~ kjij. 
v 

The Latin indices are purely spatial and k is the unit vector k/u>. The constant A 
is an energy renormalization term, having the role that in a covariant treatment is 
played by mass counter-term. 

It has been already shown [2] that the evolution operator takes the form 



U(t) = exp 



-iff D(T,T')dT,dT' A/" exp / [-iV-(r)dr], 

J o J o Jo 



(20) 
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with V(t) = e iH t Ve -iH t 

Due to the actual form of V(t), in the present case we have 



D(t,t') = Y,Di,j(t,t') (21a) 
i, j 

and 

Di,j(r,r') = <\V(V I (r)Vj(r'))\>=^J2 / ^Hl^i ' ■ v T vj ■ i Ku ■ vj 

[ e i(fc.(r / -r J ) e i(a,-fc.« J )r' e -i(^-fc.« / )r^^ r _ ^ + ^ ^ j)j 

(216) 

From the form of U one sees that the radiation of one soft particle involves the 
absolute square of the matrix element < | Vj(r) \k, v >: which is made up by addenda 
of the form: 



Ai,j(r,r') =< |VKr)|fc,i/>< k, v\Vj{t')\ > 



K 2 m 2 



r ]3h, 

" ■' - " - - --> r i(fc-(rj— rj)'\{u> — k-vj)T — i(u> — k-vi)r] 



v 



(22) 

It is only a matter of calculation to verify that 3? J o J o drdr'D = J o J o drdr'A. 

For the terms where I = J this is precisely the result already used in [2], the 
exponential of —3? J 1 J 1 drdr'D coming from eq.(20) is precisely compensated by 
the exponentiation of f* drdr'A, which comes from the sum over the radiated 
bosons, when they build up a coherent state, the only difference comes from the 
polarization factors. The case Ai :2 + A 2 ,\ involves explicitly the difference r 2 — r\ 
through the term R = k • (f 2 — ri) = uj{t 2 z 2 — r\Zi), having introduced the cosines, 
zj = cos 6>j, of the angles between k and vj. 

The actual expression we get for A = f f drdr'\A\^ 2 + A 2 ,i] is: 

A -(2^J^ dU (1 - v lZl )(l - v 2 z 2 ) 
cosujR — cosuj[R — (1 — v\Z\)t] — cosu[R + (1 — v 2 z 2 )t] + cosuj[R + v\Z\t — v 2 z 2 f^ 

(23) 

and coincides with 3? f J o dTdr'[D lj2 + D 2j {\. So here also there is a compensation 
between virtual and real contributions. 

The term coming from QD, which remains there and cannot be eliminated by 
some renormalization procedure just because it depends on f 2 — fi, is a part of the 
rescattering phase: the phase that in QED, we would call Coulomb phase *. 



* It is not, however the complete rescattering phase: since the treatment is not 
covariant there is a contribution coming from the static gravitational interaction; in 
the covariant perturbative treatment, shown in sect. II, these two terms are never 
separated. 
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We are in any case interested in the behaviour for large t, so there is a question 
on how we deal with R. The simpler attitude, followed in [2] , was to think of fj as 
independent of t, plane wave for the outgoing particles, if however we think, more 
realistically, that we have wave packets, although very broad, we would take, say, 
rj = sj + vjt, with fixed sj in the wave functions of the scattered particles, and 
so R = S + (V2Z2 — vizi)t, we see, however, that this substitution would amount in 
the previous expression eq.(23) to a substitution of vj with —vj and of R with S, 
this will not affect the t — > oo limit. 

These calculations show what really one expected i.e. that the Bloch-Nordsieck 
cancellation works also for two interacting outgoing particles; more precisely, that 
the diagonal terms compensates among themselves and the same do the cross terms. 
So now the study of the evolution of the superposition of two states of motion can 
be done, as in QED, looking only at one outgoing particle. 



B. Superposition of two states of motion 



There are two differences with the QED: one is that instead of e we find wwy, 
the second lies in the polarization term, but once we look at the evolution of one 
final particle with its radiation all the qualitative features remain the same as in 
QED, so only the initial definitions and the differences with respect to the previous 
case will be discussed. If the state after the scattering is a superposition of two 
states characterized by velocities v = and V- there are two transition amplitudes 
generated by the evolution operator U(t), say T + and 71 and then an interference 
term is produced in the transition probability at finite time: 

l(t) = 9t[T*T+] 

Each of the two factors carries its virtual correction term as discussed before and 
in [2], but the real emission of soft bosons gives, for each boson which is emitted, a 
factor: *: 

(ftm 7 ) 2 f £ z (u+-f -f-u-) 
q(ujR) = . . „ / ail — ^ =; uduj 

[l — exp[— i(k ■ v + — u)t\ — exp[i(/c • V- — u)t] + exp[i(/c ■ V- — k ■ v+)t]\ . 

/ ' (24) 

Summing over all the emitted bosons this term is exponentiated, whereas the virtual 
correction terms are noting but the negative exponentials of -^q, once with all the 
velocities equal to v+ and once with all the velocities equal to V-. 

The long time evolution produces also here a dominant term in ha(uRt) whose 
coefficient X = J dfl^ can be calculated starting from eq.(24); assuming here also 



* the real part of eq.(24) has the same form as the integrand of eq. (23), the 
quantities entering in them have however different meaning: here v + and i>_ refer 
to two different states of motion of the same particle, in eq.(23) v\ and V2 refer to 
the two final particles. 
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that the two speeds are equal and calling S the angle between the two directions it 
results: 

Gm 2 7 2 4 (cos 5 — cos0 + cos#_) 2 — |(sin# + sin#_) 2 



(1 — v cos #+)(! — v cos 9- 



(25) 



This expression could have been directly obtained from the corresponding one for 
the QED case with the appropriate substitutions of the coupling constant and of the 
spin projector, all the previous part of Sect. Ill may be interpreted as a justification 
of eq.(25). 

At this point we must extract some simpler information out of the general form. 
In the limit 6 — > one gets 



Co 



Gm 2 -f 2 4 
2tt 2 V (l-vcosfl) 5 



(sin#) z 



so that 



X. 



Gm 2 Y 2 
v 



7T Z 



2v --v 3 - (l-v' 2 )ln- 
3 K ' 1 



v J 



(26) 



This expression has, in turn, the two limits 

8Gm 2 7 2 



1 X n 



3tt 



X 



Gm 2 4 16 
v — 

7T 15 



1 + 



7 J 



For the real case 5 ^ we consider in any case the configurations at small 5, but 
with the condition $7 not small, as in [2]. After some lengthy but straightforward 
calculations we obtain 



1 X 5 



Gm 2 ^ 2 r8 
" -3 



n 



-5 2 (l-lnl(5 2 + l/ r 



The difference between X and X a expresses, as discussed in [2], the lack of compen- 
sation between the real corrections and the virtual corrections and gives the rate of 
decay of the interference term X with time *: 



Ah 



At2) 



t2 



with v = 



Gm 2 ^ 2 



n 



5 2 (l-lnl(5 2 + l/ r 



(27) 



Since we are dealing with gravitational effects we investigate the possible effect for 
large (macroscopic) bodies, in this case the speeds are certainly small, we keep 
simply the term proportional to v 4 . The resulting expression is: 



X s 



Gm 



2 2 
v 4 — (8 -7 sin 2 S) 

TV 15 



* The conditions that have been assumed for 5 and 7 make the logarithm negative 
and so v > holds. 
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And the comparison with the same expression at 6 = gives for the suppression of 
the interference term the behaviour: 



1(h) __ \h 



with v = 



Gm 2 4 




(28) 



T{t 2 ) [t 2 



TV 



IV. Some conclusions 

The results presented here are intended to be an extension of the analogous 
investigation concerning the photon emission [2]. The low-energy gravitational ef- 
fects are quite negligible * for elementary particles, anyhow the explicit computation 
shows a time decay of the interference term which is exponential in t, the factor 
that multiplies the time is different in the two cases: its energy dependence is much 
stronger, the angular dependence shows some minor differences, if fact both in grav- 
ity and in QED the angular dependence is logarithmic for small angles. It is however 
quite plausible that the same mechanism is at work for macroscopic bodies where 
it could became comparable and ever larger than the electromagnetic effects. The 
actual expressions intended to refer to macroscopic bodies, eq.(28) may be ques- 
tioned on the basis that the starting point was a Bloch-Nordsieck model designed 
for high energy collision, looking in detail to the procedure that has been used one 
sees that the essential point of the approximation was to neglect the change of ve- 
locity of the emitting body in the process of radiation, now it is clear that for a 
long wave-length emission having a massive body is enough to make the variation 
of the velocity negligible also for small velocities. In practice interactions with the 
rest of the world will be, for a macroscopic body, certainly more relevant than grav- 
itational radiation, but this kind of radiative process last should be quite universal 
and for this fact it has some interest of principle. Another point to be stressed is 
the possibility that this kind of dynamics acts also on photons. 

As it has discussed in [2] this particular decoherence effect is evidently due 
to the presence of massless particles, and so to the possibility of producing in the 
collision process of an indefinite number of them, so every exclusive channel gets, 
at the end, probability zero; in an S-matrix treatment the sum over different final 
states is not simply a practical procedure, but also a necessity in principle. As the 
radiated particles become softer and softer the time needed to radiate them becomes 
longer and longer, for every finite time no true i.r. divergence is present [9], the 
zero value for of the probability for every exclusive channel in anyhow attained with 
continuity in time. 

There is a more complicated effect of decoherence, which affects in case of 
superposition of two (or more) states of motion the rescattering phase, see eq.(7d), 
but the possibility of observing these effects appear even more remote and also less 
relevant in principle, so it has not been computed. 

It can also be noticed that the result does not comes purely from kinematics, 
the particular form of the coupling is relevant, if e.g. one would consider a chiral 

* needless to say: the high-energy gravitational effects are unknown 
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coupling, the insertion of a massless pseudoscalar particle on an external line would 
give 

1 k 

-75«(p) = ^ i pjIMp) 



p-k-m 2p-k- k 1 

and there would be no divergence in the limit k — > 0. 

Tentatively one could argue that the emission both of photons and of gravitons 
arises from a dynamics giving a privileged role to the momentum variables and 
therefore the persistence of states which are superposition of two or more sharply 
different momenta is not tolerated. We know also that this kind of emission is 
related to the presence of long-range forces [3], it is not obvious whether this has 
to do with the decoherence process, but it is certainly the reason why both the 
electromagnetic and the gravitational interaction sum up starting from microscopic 
size to macroscopic size and thus making a qualitative continuity in the dynamics 
of small and large objects. 

The effects of the infrared radiations are clearly due to a dynamics which is 
quite standard; this can be said beyond any doubt for the electromagnetic infrared 
emission, which is well confirmed experimentally, but also for the gravitational 
radiation, for which a direct experimental verification is not available, the common 
lore is to take for granted its existence. So these effects should be there even if the 
quantum mechanics should undergo some important modifications which, in very 
different forms have been postulated also in the aim of introducing a decoherence 
while keeping the relevant features of quantum physics. This is not the place to 
discuss the very large literature on this field, as an example one of the most carefully 
elaborated proposal is mentioned [10]. There a modification in the in the evolution 
of every quantum system is postulated, by introducing a non-Hamiltonian term, 
while the Hamiltonian part is not modified; in the present treatment the evolution 
remains Hamiltonian, but in the limit t — > +oo the evolution operator is no longer 
unitary, while remaining isometric. 
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Figure captions 

1. ) The basic 2^2 scattering graph, for particles p,p' going into particles 

q,q'- 

2. ) The four graphs for 2^2 scattering with emission of a soft particle with 
very low four- momentum k. 

3. ) Three of the six graphs for 2 — > 2 scattering with virtual correction in s, t, u 
channels, the virtual particle bears the four-momentum /. 

4. ) One of the four 2^2 scattering graph with self-energy correction. 
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